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TECHNICAL NOTE 2^73 


ON THE SPECOEIM OF ISOTROPIC TURBULENCE 
By H. W. Llepmann , J. Laufer, and* Kate Llepmann 

SIWMARY 


Measurements of the spectrim and correlation functions at large 

Reynolds number ^RN » lo5 "based on the grid mesh^ have been made, as 
well as a series of accurate spectrum measurements at lower Reynolds 

number ^RN « 10^). 

The results are compared with the theoretical laws proposed In 
recent years. It Is foimd that the measurements at large Reynolds 
numbers exhibit a range of frequencies where the spectrum is nearly of 

the form n”^^^ , 

The largest part of the spectrum in the initial stage of decay at 
the lower Reynolds number was found to follow closely the simple spec- 
trum — — , where A and B are constants and n is the frequency 

B + n2 

of fluctuation. At = 1000 (where x is the distance behind the 

the grid and M is the mesh size) the spectrum approaches a Gaussian 
distribution. 

The second, fourth, and sixth moments of the spectrum have been 
compu"ted from the measurements euad are discussed in relation "to theo- 
retical results. 

The significance of the number of zeros of the fluctuating velocity 
u(t) is discussed and examples of measurements for the determination of 
the microscale of turbulence X from zero counts are given. 

INTRODUCTION 


A field of turbulent fluctviations represents a dissipative system. 
Viscosity will tend to dissipate the energy of the fluctuations into 
heat and, if a stationary state of the system is to be kept, energy 
must be continuously supplied to the fluctuations. It is characteristic 
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of fluid flow phenomena that the rate of dissipation is strongly influ- 
enced hy the inertia terms in the equation of motion. These convective, 
nonlinear terms do not describe any mechanism for the production of heat 
from the kinetic energy of the fluctuating field and their influcence 
upon the rate of dissipation is therefore indirect. The rate of dissi- 
pation is proportional to the mean square of the vorticity. The non- 
linear inertia of the fluid motion usually tends to increase this mean 
square. A similar mechanism is well-known in fluid mechanics in the 
shock-wave formation where the tendency to steepen the shock front is 
balanced by the tendency to diffuse it because of the viscosity and heat 
conduction of the fluid. 

Ohe understanding and smalytlcal formulation of the nonllneeu: effects 
upon the vorticity distribution is a central problem in turbulence 
research. In the search for an understanding of this mechanism Taylor’s 
Introduction of the concept of isotropic turbulence was very Important 
(reference l) . Isotropic turbulence represents a much s inkier type of 
turbulent motion than the general shear flow problem but does Include 
many characteristic features of the general problem. Most of the recent 
progress in turbulence research came from investigations of Isotropic 
turbulence. 'K£rm£xL (reference 2) introduced the convenient” formalism 
of the correlation tensors of different-rank. KtCrman and Hbwarth 
(reference 3) then gave the first equation in a suitable form relating 
the double correlation function to the triple correlation function. 

This equation shows the balance between inertia and viscous forces in a 
turbulent field. The triple correlation represents the nonlinear terms 
of the equation of motion. The K^rman-Hbwarth equation is clearly 
indeterminate as long as one does not have an additional analytical 
representation, in terms of the double correlation functions, of this 
fundamental process of vortex regrouping. 

The theoretical research of the last few years has made some pro- 
gress at least toward an understanding of Isotropic turbulence. 

Kolmogoroff (reference 4) introduced the concept of local isotropy into 
turbulence, that is, the hypothesis that the motion of the smaller eddies 
in turbulent flow is always isotropic. This hypothesis, for which there 
is some experimental verification at present, makes an investigation of 
isotropic turbulence even more interesting. Kolmogoroff then proceeded, 
essentially on the basis of dimensional analysis and very simple physical 
reasoning, to arrive at some results concerning the correlation function 
and so forth. A thorough review of Kolmogoroff *s work has been given 
by Batchelor (reference 5)* CJnsager (reference 6), Weizsficker (refer- 
ence 7), and especially Bielsenberg (reference 6) attacked the problem from 
the point of view of the turbulent spectrum rather than the correlation 
function. Taylor had Introduced the concept of the spectrum into 
turbulence as early as 1938 (reference 9). However, the one- 
dimensional form of the spectrum as given by Taylor is very convenient 
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for experimental measurements tut not too convenient for theoretical 
study. Heisenberg therefore lntrod.uced. a three-dimensional space spec- 
trum function E(k) . The function E(k) dk Is defined as the fraction 
of the total turbulent energy with wave number between k and k + dk. 
Taylor’s spectrum S^(k]_j dk^^ Is defined as the fraction of turbulent 

energy with components of the wave-number vector k in one fixed direc- 
tion between kn and k-i + dk, . Heisenberg's and Taylor's spectrum 

functions are Fourier transforms of the correlation functions of Karman 
Corresponding to the K^umian-Hbwartb equation, a similar equation for 
E(k) may be written which obviously does not yield anything different 
from the equations for the correlation function. ■ However, It Is possible 
that the concept of spectzrum Is somewhat more Intuitive than that of 
correlation function - probably because spectrum Is a more familiar 
concept from other fields of physics - and most recent work operates 
essentially with the concept of spectrum. Knowledge of the spectrum 
function E(k) Is sufficient to describe the turbulent field ccMnpletely. 
A rather complete survey of this group of recent theoretical Investiga- 
tions and the experimental evidence can be found In Batchelor's lecture 
to the Seventh Intematlonal Congress for Applied Mechanics (refer- 
ence lO) . A brief discussion of some theoretical results for later 
reference Is included In the section "General Considerations" of this 
report . 

The general aim of the present experimental Investigation Is a 
study of the nonlinear exchange mechanism. There Is not at present any 
straightforward and simple way to measure the exchange terms directly. 

In fact, the research has to start on a more basic level. A fairly 
broad Investigation Is needed to make sure that the flow conditions are 
close to those assumed by theory and that the experimental equipment Is 
satisfactory. Finally, a rather thorough knowledge of various aspects 
of the turbulent field Is needed before experimental results can be 
evaluated Intelligently. 

In the present work the following points were Investigated: 

(1) The correlation function and the spectrum at large Reynolds 

numbers 

(2) The decay, microscale, and spectrum at Intermediate and low 

Reynolds mmibers 

(3) ^Ehe zeros of the fluctuating velocity components and their 

relation to the spectrum and moments of the spectrum 

Fart of these measurements can serve as an Independent check of 
results obtained by Batchelor and Townsend and others, and part of the 
results are new. It should be emphasized here that the ultimate aim of 
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a more or less direct- measurement of the exchange term has not“been 
carried out”^t and, in factT^was not even attempted. It was felt neces- 
sary first to improve and refine the methods’, especially the measurement 
of the s^ctrum, md, as mentioned before, to gain a view of many aspects 
of the turbulent field. 

The inveatigations were carried out at the Guggenheim Aeronautical 
Laboratory, California Institute of Technology, as part of the turbulence 
research conducted under the sponsorship and with the financial assist- 
ance of the national Advisory Committee for Aeronautics. The authors 
would like to acknowledge the cooperation of Mr. F. K. Chuang Eind the 
stimulating discussions with Professors Lager strom and DePrlma. Mr. M. 
Jessey was responsible for the design and construction of most of the 
electronic equipment and this essential help is gratefully acknowledged. 
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SYMBOLS 

turbulent velocity fluctuation (1=1, 2, 3) 
space coordinate (k = 1, 2, 3) 
mean velocity 

mean square of velocity components 
Karman's double correlation tensor 
Batchelor’s spectral tensor. 

wave -number vector 
components of" wave -number vector 
frequency of fluctuation 
Heisenberg's spectrum density 
Taylor’s spectrum density in space 
Taylor’s spectrum density in time 
dljienslonleBs spectrum density parameter 

dimensionless frequency parameter 
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Ry = f(r) 


double correlation, functions 
triple correlation functions 
scale of turbulence 
microscale of turbulence 

average number of "I values” per unit time of a random 
function 

viscosity 

klnanatlc viscosity 
dimensionless parameter 

= r distance between two points In field 
value of random variable u(t) 
value of random variable du/dt 
density 

pressure fluctuation 
spherical coordinates In k-space 

for case where r Is measured In direction norma], to Uj^ 


P(N) 


rate of turbulent energy dissipation 
probability distribution of N| 

GSnOEAL CONSIDEBATIONS 


Definition of Terms 


Homogeneity, Isotropy, and stationary state .- A field of turbulent 
fluctuations in space Is given by the velocity u = j^u^^Xyt^J . The 
field Is called homogeneous If all mean values of u anfl their 


V 
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derivatives do not depend on the value of Xj. It Is called Isotropic p 

If these mean values are Independent of reflection and rotation of- the 
coordinate syst«n. In most measurements of-^;urhulence a velocity com- 
ponent Uj^ Is observed as a function of time at a fixed position. For ^ 

example, measurements are made In a vlnd tunnel at some distance behind 
the grid. In this case the process Is called stationary^ mean values 
of the square of u(t) emd Its derivatives do not depend upon time and 
the correlation functions depend only upon the relative time Interval. 

In the definition here one has to use the ensemble average to take' a 
mean value. If the ordinary time average Is used, the statement that 
the mean values do not depend on time becomes trivial. A stationary 
process In time Is analogous In definition to homogeneity In space. 

Isotropy has no such counterpart since time is a single coordinate. 

Time and space derivatives .- It Is typical of present turbulence 
research that measurements are often made In a moving stream at a fixed 
position. On the other hand, theoretical considerations often deal vlth 
space distribution of a turbulent field at-rest. To relate experiments 
and theory one uses the fact that the space coordinate In the direction 
of the mean motion Is tlmellke. This fact is sometimes e:qxres8ed by 
stating that the operator &/bt can be replaced by -U hfdx if U Is 
the mean velocity In the direction x. t** 

This statement is obviously never correct rigorovisly. In certain » 

applications the difference, that— is, the error, cein be given or esti- 
mated. For example, for any stationary process u(t) which possesses 
a derivative 


“ i - ° (1) 

applying the exchange of operators simply gives 

u = -uTg = -U 

But 

-U ^ 0 

OK 

since this term determines the decay of turbulence and for isotropic 
turbulence has the value 

-U ^ = lOv^ 

Sx ^2 T 
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It Is possible, however, to use the Interchange between time and space 
variables In certain mean values. For exan^le. It has been found experi- 
mentally by Townsend (reference U), In agreement with the measurements 
presented here, that the mean square of the time derivative agrees 
closely with the mean square of the space derivative In the direction 
of meGUi motion; that Is, 



( 2 ) 


To estimate here the degree of approximation within which this 
result Is tzrue Is more complicated. One can write down the equations 
of motion and compute these terms and then try to estimate the differ- 
ence; for example, for 


Hence, 


with 



P OXi 



+ B 




The largest term In R appears to be 




du4 du4 




\2 


Heglectlng this term compared with 


4 ?^ 

W/ 


error of the order of the turbulence level 


2 

I hence would amount 



to an 


Mean values .- So far mean values have been used without any defini- 
tion of what Is meant by putting a bar over certain quantities. As a 
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matter of fact, the definition and use of mean values In turbulence Is 
not an easy problem. Experimentally, the average Is, In general, the 
time average provided by the response of a measuring Instrument, for 
example, a thermo-cross-galvanometer combination. OSils time average Is 


equivalent to 



0 dt '(diere t Is time and 9 


Is the time interval 


for the averaging. For theoretical considerations an ensemble average 
is often more convenient. ExpertmentfOly^ an ensemble average can be 
obtained by averaging over single measurements separated by time Inter- 
vals STofflclently large to have zero time correlation. If the time 9 
In the averaging process is large compared with the correlation time 
scale T, then ore should expect the time and ensemble averages to agree. 


Ihe finite extent of the measvirlng probe, that is, specifically, the 
length of the hot^vlre. Introduces an additional average over space In 
turbulence measurements. If the length of the hotiawlre Is not too great 
this can be accounted for by making corrections based on a rou£^ estimate 
of the space characteristics of the turbulent field. 


Correlation functions, spectrum .- The space distribution of a 
homogeneous turbulent field is characterized by giving the space correla- 
tions in the form of the well-known K^rm^ correlation tensors. 





r = 






>- 


(3) 


The space spectrum is then most conveniently defined as the Fourier 

transform of (references 12 emd I 3 ) . Let k denote the wave- 

number vector, that is, the vector in the direction of wave propagation 
with absolute magnitude k = |k| = 2 x/a, A being the wave length. 

Then a spectral tensor cein be defined by 




Kl3(r)e^(^-S) dr(r) 


(4) 
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Heisenberg's spectrum E(k), that Is, the energy density of turbulent 
fluctuations contained between k and k + dk, is then obtained as 
the mean value 


E(k) 


^ pt p2« 

Jo Jo 


sin 0 d0 dcp 


(5) 


For isotropic turbulence and ate each expressible by a 

single function f(r) and E(k), respectively. 

The spectrum function P(n) Introduced by Taylor was originally 
defined in time frequency n rather than in space. The relation 
between F(n) and the Heisenberg spectrum thus involves an interchange 
between time and space variables. The equivalent definition of the 
Taylor spectrim in space F(kQ_) is the fraction of turbulent energy 

contained between kj^ and ki|^ + dk^ where k^ denotes the component 

of k in a fixed direction. The direction is, in general^ chosen to 
correspond to the direction of mean speed U in the experimental setup. 

There exists a simple relation between F, the Taylor spectrum 
function, and E, Heisenberg's spectrum function, wbicb is valid for 
homogeneous, isotropic turbulence. 




E(k) dk = 1 


( 6 ) 


r 

w 


where from equation (5) 


P(ki) = ^ f(r)cos k^r dr 


f(r) “j^ P(k 2 ^cos k^r dk^ 




(7) 
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If one considers again the stationary time process, then one can define 
the corresponding time correlation and time spectrum as Is standard 
practice. If u(t) Is a certain velocity component measured at one 
space position In the course of time, then one can define 


^(t) ^ (8) 

u2 

where t denotes a time interval. The corresponding spectrum F(n) 

is then defined as the fraction of u^ contained in the frequency hand 
n, n + dn, that is, the so-called "power spectrum." One has th^: 


J pOO 

f i(T) COB 23 tnr dr 
0 


poo 

i|f(T) = I F(n) cos 2jtnT dn 

From equations (7) and (lO), respectively, follow then the relations 
between spectrum, correlation function, and mean derivatives of the 
fluctuation which are so often used in turbulence. Specifically, the 
following qviantities will be used: 

(l) The space scale L and time scale T of turbulence defined by: 


(9) 


( 10 ) 


I 


L 



dr 



T = dr 


F(0) 

k 


( 11 ) 


r- 

w 


( 12 ) 
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(2) The so-called microscale X or tte corresponding mean-square 
frequency n2 defined by 




dk^ 


(13) 


^ n%(n) dn 

From equations (7) and (lO), clearly: 

= «f«(o) 

■m 




(15) 


4A2 = -t"(0) 

■iSf 


( 16 ) 


Triple correlations .- For a description of the dynamics of turbu- 
lence one needs also the triple correlation functions vhich are intimately 
connected with the nonlinear terms in the equations of motion. For iso- 
tropic, homogeneous tvirbulence the triple correlation tensor of Karmen 
can be expressed again by a single function, for example, h(r) vhere 


(^) 3/2 


Cl’ =x^ + 6^ 


’ = X. 


1 = 2,3 


(IT) 
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Functions equivalent to the spectral densities can he Introduced 
here as vas flrsh done by Lin (reference 13) • For example, one can 
define a function equivalent to the Taylor spectrum for the 

double correlations and vrlte 





sin k^r dr 


(18) 


Similar relations In time can also be -written. 


Equations of Motion 

The equation for the correlation function In homogeneous. Isotropic 
turbulence Is the Keumilln-Howarth equation which Is, of course, a conse- 
quence of the Havler -Stokes equa"tlons: 

^ . i |£) (19) 

St ^ ' \Sr r/ rSrj 

where f and h are the double and triple correlation functions defined 
abo-ve. 

The corresponding equation In the -wa-ve-number space has been given 
by Lin (reference l4) ; 


+ W(k,t) = -2vk2E(k) 


( 20 ) 


with 

W(k,t) = |t2jk2Hi'*(k) - kHi'(k^ 

A thorough discussion of these equations Is given by Karman and Lin 
(reference 15) . 


Form of Spectrum Functions 

The aim of a theory of Isotropic turbulence Is the determlna-tlon of 
the function f(r) or E(k) . Earlier attenqpts to ob-taln f(r) are 
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4 found, for exaDQ>le, in the Kerman and Howarth paper (reference 3)* In 

recent years it has hecome more common to operate with the spectrum 
function rather than with the correlation function. 

u 

In the present state of theoretical analysis of E(k) or f(r) 
it is best to distinguish tvro steps. First, one can arrive at some 
conclusions regarding the general form of E or f in certain regions 
without making any explicit statement about the nonlinear term W(k,t) 
or h(r) . Second, one can arrive at conclusions by assuming a specific 
form of W or h. 

The first set of theoretical results includes Kolmogoroff’s theory 
of locally isotropic turbulence and the results obtained concerning E 

_5/o 

and f, especially the k ' law for E(k) . This has also been 
obtained Independently by Welzsacker (reference 7) and (kisager (refer- 
ence 6) . The set Includes also the k^ law for small wave numbers 
obtained by Lin (reference 1^), Batchelor (reference 12), and others 
and, further, Loitslanskll's (reference l6) result concerning the invari- 
ance of u2 / r^f(r) dr during decay. Finally, the behavior of the 
I Jo 

turbulent field at very low Reynolds numbers as already discussed by 
ii K^irm^n euid Howarth, Loltalanskli, and recently again by Batchelor and 

Townsend (reference 17) belongs to this group. 

On the other hand, the k law of Heisenberg (reference 8) for 
leurge wave numbers and all relations obtained for medium wave numbers 
depend, on assumed forms of the exchange term W(k,t), h(r,t), or 

corresponding terms. K£rm&i (reference l8) has recently systematized 
these assumptions . An analysis of all these various laws and assumptions 
can be found in Batchelor's lecture to the Seventh International Congress 
for Applied Mechanics (reference 10) and in the recent paper by Tra.mtfm 
and Lin (reference 15) . 

From the point of view of experiment these two goups of results 
differ essentially. Verification of results of the first group mainly 
means a verification of rather broad and essentially kinematic assump- 
tions . For example, the theory has until now always been based on the 
consideration of an unlimited turbulent field with zero mean velocity. 
Experiments are carried out with turbulence which is produced by a 
specific mechanism and which is certainly somewhat inhomogeneous turbu- 
lence. Furthermore, the nwvlnnnn length or mlnlmiTni wave number in an 
experimental setup is limited. The range or validity In Reynolds number, 

for example, for the k“^/^ law, is not easily predictable and, similarly, 
the range of wave numbers for the k^ law. 


<a 
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Consequently, the aim of experiments here Is to estahllsh the 
applicability of the theory. Formally, there is little doubt that these 
laws apply in the assumed klnematical model since the derivation from 
the equations does not Involve any additional assumptions. 

The aim of experiments for the verification of the second group of 
results goes deeper because, in addition to the applicability of- the 
theories to the experimental setup, the physical mecheuilsm of exchange 
is involved here. It is consequently impossible to work on the second 
group of problems fruitfully without having gone to a certain extent 
through the first. 

A further question to the experimenter is the consistency of measure- 
ments. Quite often the theoretically defined quantity cannot be measured 
directly. Finally, the sensitivity of the form of- the measured quantity, 
as, for example, l(k), toward the assumptions meuie is Important. H3ils 
often makes a decision In favor of a certain specific assimqjtion very 
dubious . 

These simple considerations are included here to Justify the 
approach to the problem. Even if the ultimate aim is quite clearly 
recognized it is not always possible to direct" the experiments Immediately 
toward this aim since any result- will be ambiguous. Clearly, a direct 
experimental approach to the turbulent exchange problem is most desired 
but this approach does not seem very fruitful before the background is 
somewhat cleared up. 


Remarks on Turbulent Fluctuations as Stationary Random Processes 

If a turbulent velocity component u is measured at a fixed posi- 
tion as function of- time, u(t) denotes a stationary random process. 
Processes of this general type have been considered in physics and math- 
ematics before and it is InteTOstlng to see how turbulence fits into 
these known processes. Brownian motion and the shot effect are examples 
of processes of this type. Often the term "random noise" is applied to 
processes of this kind. If one measures the turbulence spectrum in the 
standard fashion by analyzing the output of a single stationary hot-wire, 
or if one measures the time derivative in the same way by differentiating 
the output current of a single hot-wire, and so forth, in all these cases 
one analyzes a statlonaxy one -dimensional random process and application 
of these results to the three-dimensional space field of turbulence 
requires additional steps, such as the Interchange of space and time 
variables discussed above. 

Before the results of measurements are given in this report it 
appears useful to discuss briefly two aspects of the problem of 
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stationary random process which are used in the measurements, namely, the 
power spectrum and the number of zeros of u(t) per unit time. 

Power spectrum .- If u(t) denotes one turbulent velocity component, 
not necessarily the one in the direction of mean motion, one can define 
a power spectrum F(n) as discussed in the section "Definition of Terms." 
The necessary conditions underlying the concept can be found in the 
literature, for example, in references 19 and 20. 

Dryden (reference 2l), who cairied out one of "Hie earliest measure- 
ments of the turbulent spectrum, found that the measurements of F(n) 

• could be represented rather well by the function 


F(n) 


F(0) 


1 + c^n^ 


( 21 ) 


where c and F(0) are characteristic constants related by the normaliza- 
tion of F(n). The corresponding time correlation i|r(T) la then 


+(t) 



( 22 ) 


Clearly, a function of this form cannot represent the spectrum for 
arbitrary high values of n since then 


n?F(n) dn = 00 


Thus the dissipation of turbulence - if for the moment one takes the 


for granted - becomes infinite. Indeed, Dryden *s 

measurements showed that for large values of n the measured spectrum 
falls below the values given by equation (21) . 

The measurements of F(n) made during the course of the present 
Investigations are accurate enough to obtain 




n 


2 



n%(n) dn 


the second moment of F 



l6 


NACA TN 2l^73 


as veil as 


pQO 

^ = / n^(n) dn 

Jo 


the fourth moment of F 


both of vhlch are found finite - as expected. - and. In general agre^nent 
with measurements of~Batchelor and ^Downsend of the mean square of the 
first and second derivatives of u(t) . HovSver, in agreement vith 
Dryden it vas found that by far the largest part of the turbulent energy 
is contained in a portion of the spectrum F(n) for which the simple 
form of equation (21) gives a good approximation. Hence, one may take 
at present this form of F(n) as, at least, an empirical interpolation 
formula from experiments and one may Inquire into the general nature of 
such a spectrum. Without implying too much about the physical signifi- 
cance of the simple spectrum function the following properties are 
interesting and possibly significant. Dryden noted that this spectrum 
plays a specific role in random processes. Today one can formulate the 
significance of this spectrum as follows: 

(1) If the probability distribution of u(t) is Gaussian - as is 
known with good accuracy from experiments of Simmons and Salter (refer- 
ence 19) and Townsend (reference 11) - and the spectrunrhas the form 

— ^ , then the process is Markoff ian; that is, it is the simplest 

B + n2 ■ ■ ... 

random process for a continuous variable (see reference 20). 

(2) Processes leading to a spectrum of this type are well-known and 
fairly common. For example, the so-called random telegraph signal, that 
is, a variable vhlch attains the values a and -a Jumping from the 
one value to the other at random times, has this spectrum. The number 
of changes of sign is then given by a Poisson distribution. 


Hence, there may be a possibility of relating some part of the 
turbulence to the random shedding of vortices by the grid. This appears 
worth a further investigation. Furthermore, it is interesting to note 
that - assuming again that the exchange of_tlme and space variables is 
closely valid - the simple spectrum satisfies the kinematic requirement 
for low wave numbers. iHie corresponding Heisenberg spectrum behaves 

like k^ for k — >0. Indeed, it follows simply from equations (6) that 

the Heisenberg spectrum E(k) correspondli^ to F(n) has 

1 + c n 


the form 


E(k) 


(1 + L^k^jS 
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where L is the scale of turhulence and hence E(k) = as k — >0. 

Furthermore, for large -values of k, that Is, kL » 1 



L 

which makes it experimen-tally very difficult to distinguish this spectrum 
behavior ftom the k~5/3 law. 

It is not Intended here to overemphasize the impor-tance or physical 
significance of the single spectrum correlation functions. It is only 
pointed out that the simple spectrum is an excellent Interpolation 
formula with possible physical significance which has to be studied 
further and that one should be quite careful in the experlmen-tal verifi- 
cation of power laws such as the k“5/3 law. Using the simple spectrum 
Loltslanskli’s in-yarlant becomes simply proportional to 

u^[ot( 0)] 5 as ^(ul)5 (23) 

a form which can be experimen-tally determined. Results are gl-v-en later 
in this report. 

Number of zeros of u(t) .- Another rather interesting problem 
arising in s-tationary random processes is the expected number of passages 
through zero (or any other cons-tant -value) per unit time. Questions of 
this nature have been discussed in recent years, for example, by Rice 

(reference 22).^ 

Let P(l^Ti) d| dq denote the probability at a given time of finding 
u(t) be-tween the -values 5 and | + d| and du/dt be-bween q 
and Ti + dq. Then p(S,ii) d§ dq denotes the time the trace spends in 
the Interval | + dg having a derlvatl-ve between ti and tj + dq. 

To ob-tain the number of passages -through the interval dg, regardless 
of the time spent at a crossing, p(g,T]) d| dT] has to be divided by the 



^Turther references -to other investigations In-to the subject, espe- 
cially the work of Eac, are given in Rice’s paper. 
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time spent in the intisrval at a crossing to reduce all passages to the 
same statistical weight. Let t he this time; then cleeo:ly 

d| 

"'In 

since is the "velocity" in crossing the . interval d| and | t] | has 
to be used since t does not depend on the sign of t| . Hence , the 
number of passages through | + d| with a specific value of the 

derivative tj Is given by 

and, consequently, the total number of passeiges through the interval 
regardless of ti, by 


N| |n|p(&>Ti) dri (2l^) 

This is the formula given by Rice in reference 22 where a moire rigorous 
proof wllL be found. 

So far the assumptions underlying N| are very general. To obtain 
a simple relation between the spectrum and additional assumptions 

are necessary. These assvonptions are jiistlfled for linear processes, 
such as electrical noise problems and so forth, but certainly not 
rigorously time for tvirbulence. Still it is of- considerable Interest to 
study the effect of possible deviation from these assumptions on N|. 
Since u(t) is a stationary random process 


du 


u 


dt 


0 


that is, there is no correlation between u and du/dt. If both u 
and du/dt have Gaussian distributions, then it is directly follows that 
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Specifically, the mmber of zeros Nq 1b given hy 



(26) 


In this form formula Nq quite similar to the harmonic value for 
which the number of zeros Nq Is given hy 


Nq - ^ 


The distribution of du/dt is, however. In geiieral, not Gaussian, This 
fact is essentially a consequence of the nonlinearity of the equations of 
motion and was verified e^qjerlmentally hy Townsend (reference ll) . Still, 
It is interesting to Investigate the deviation of Nq frcan the mean- 
square frequency as obtained hy other means, for example, from P(n) 



and to extend this comparison to Reynolds numbers so low 


that the Influence of the nonlinear terms should become unimportant. 
Unfortunately, an acciirate comparison cannot he given as yet. It turned 
out that It is experimentally easy to count the number of zeros with 
fair absolute accuracy, that Is, of the order of ±10 percent. To 
Increase the accuracy beyond this point proved to he more difficult 
because of the very large range of passage time which the counting 
apparatus had to reproduce accurately. Consequently, only preliminary 
measiirementB are briefly reported here. 


APPARATUS AND METHODS 
Wind Tunnel' 


The wind tunnel used for most measurements was the small 20- inch- 
square GALCIT "Correlation" tunnel built many years ago for similar 
investigations. Repeated Improvements in the tunnel have reduced the 
free- stream turbulence level to the order of 0.03 percent. Traversing 
equipment along the tTjnnel axis for about 200 centimeters' length is 
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provided. The GALCIT 10-foot tunnelLvafl used for some investigations 
at high Reynolds numbers reported here. Slcetches of the tvo tunnels 
are Included as figures 1 and 2. 


Grids 

The grids vised in the present investigation differed somewhat from 
previous ones. They consisted of precision^ voven-hrass screen mounted 
in wooden frames. The screens, obtained from Tyler and Company in 
Cleveland, are very uniform. The geometry of the grids is given below. 


Grid 

Rod diameter, d 
(cm) 

Mesh size, M 
(cm) 

1 

0.4i3 

1.68 

2 

.184 

.818 

3 

.108 

.422 

4- 

.044-5 

.l4l 


Only in the 10-foot tunnel was a grid of the two-plane type used. The 
dimensions for this grid were M ■ 4 inches and d » 0.T5 inch. 


Hot-Wire Apparatus 

Except- for the investigations in the 10- foot tunnel, wires of 
0.0000^- inch diameter were used. In the 10- foot tunnel the wire diameter 
was 0.0002^ inch. The wires varied in length, depending on the problem, 
between about 0.3 centimeter and 0.1 centimeter. The silver cover of the 
platinum wire was in all cases etched off before soMerlng the wires to 
the tips of sewing needles. 

The amplifier used was rebuilt for the investigation. The frequency 
response is flat within a few percent between l/2 cycle and 25,000 cycles. 
A cut-off filter at approximately 10^000 cycles is provided to reduce 
noise and pickup whenever there is no need for an investigation of higher 
frequencies. This was the case most of the time in the present set of 
experiments . 

A resistance-capacitance compensator and a stage of differentiation 
using a second compensated amplifier is incorporated. The output is 
read with a thermocross potentiometer arrangement. The compensation is 
adjusted by means of a square wave arrangement. 
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« Frequency Analyzer 

A wave analyzer from Hewlett and Packard wltli an adjoistatle "band 
- width from nominal 30 cycles to 1 U 5 cycles was used for spectrum measure- 

ments. The analyzer was slightly modified to allow reading the output 
on the thermocross meter. 

The measTorements In the 10- foot tunnel were made with a General 
Radio Sound Analyzer which is not so sMtahle for turbulence measure- 
ments as the wave analyzer. Part of the very marked Improvement in 
accuracy between earlier spectrum meas'urements in the 10- foot tunnel 
and the recent meastiremente in the 20- inch tunnel is due to Improvement 
in the equipment and part is dvie to the fact that the flow in the smaller 
tunnel is steadier than that in the 10-foot tunnel. 


Zero Counting 

The number of zeros or | values was first counted with the help 
of a photomultiplier cell. The cell was mounted behind a narrow silt in 
front of an oscilloscope screen and the passage of the trace on the 
t screen thus produced pulses which, after suitable amplification, were 

counted by means of a 2^ scaler counter. Every 512th pulse was thus 
^ registered on a mechanical counter. An average count comprised about 

1(P passages. Recently, an electronic g^te circuit was ueed for the 
same purpose. 


P 

Length Correction 

p 

The measured values of were corrected for the length I or 

the hot-wire by means of the relation ■ ( ^^)nieasured “ ^ 

where G = X^f^^(O) . The values for G given in this report are 
uncorrected since their limited accuracy does not warrant an elaborate 
correction involving still less known higher derivatives of f. 

A len^h correction should be applied to the spectrum measiirements . 
However, it appears difficult to obtain an analytical expression for such 
a correction. Check measurements were therefore made with wires of 


^Length correction formulas were derived by the present authors in 
a previous unpublished report submitted to the MCA. The procedure is 
straightforward and similar derivations have been published in the mean- 
time by Frenklel. Thus there is no need to give the derivation here. 
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different length and the correction was evaluated semlempirlcally. It~ 
turned out that the correction is unimportant for F(n) Itself for 
the wire length used. For the higher moments the correction becomes 
more Inqjortant but at the same time the accuracy of the measurements 
decreases. Thus, it was felt- that, correction should be axqplled at 

present only to The rest of the data are presented uncorrected for 

the length of- the wire. 


RESULTS OF MEASUREMENTS 
Decay 


A few measuranents of the decay of turbulence behind a grid were 
made and figure 3 presents a sample. The decay measurements agree with 
Batchelor and Townsend's (reference 23) measurements in the Inlt^l 

stage of-decayj they show the nearly lineeo: relation between l/u^ and 

t or x/M for —< 100 (approx.) and the subsequent Increase in the 
M 

rate of decay. 


Measurements 

The value of during decay was determined by Townsend's method 
from the mean square of the differentiated wire output mating use of the 
interchange between time and space derivatives. Then 


1 

^ ■ ;?W 


w 



Also was determined from the measured spectrums. The relation here 
is 



O 

' ^ Jo 
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f A set of results Is shown In figures 4 and The estimated errors are 

Indicated by the extent of the lines from the points. The straight line 
with slope lOv/u Is plotted In the figures to allow a comparison with 

* the measurements. It Is seen that the measured dX^/dx Is near lOv/U. 

The measured slope Is a little smaller than lOv/u eind the agreement Is 
not so close as that given by Batchelor and Townsend. However, In the 

Initial stage of decay there Is no doubt that dX^/dx la much closer 
to lOv/U than to 7vA^. The slope lOv/u Is that which follows from 

the assumption of Isotropic turbulence If the relation of l/u^ Is 
linear with x/m or t. The slope 7v/u follows from certain, other 
similarity assumptions (e.g., reference 24). The question of the slope 

of X^(x) has been discussed sufficiently by Batchelor and Townsend 
(references 23 and 25), Frenklel (reference 24), and others, and there 
Is no need to repeat the argrmients here. 


Spectrum and Correlation Measurements In GALCIT 10-Foot Tunnel 

More extensive measurements were made of the spectrum of turbulence. 
These measurements Include earlier Investigations of the spectrum at large 
Reynolds ntimbers In the 10-foot tunnel. These latter measurements were 
done In order to check on the existence of a range of frequencies for 
-5/3 

which the k spectral law applies at Reynolds numbers which can be 
reached In a wind tunnel. As mentioned before, these measurements do 
not have the precision of the later spectrum measurements at lower 
Reynolds numbers. The resxzlts are shown In figure 6 plotted In a reduced 
scale 

1^=0 against |F(o)n = ^ 


The simple spectium 


a = -A- 
1 + ^2 

Is shown for con 5 >arlson. The scatter of the points Is evident and any 
conclusions therefore are not too convincing. Here the simple spectnmi 
does not agree so well as in the measurements at smaller Reynolds 
numbers. Figure 7 shows the same measurements In a logarithmic plot. 

-5 /a 

The line with a slope corresponding to n Is Included. These same 
measurements were used by Kirm^n (reference 22) for a similar comparison 

with the n"^/3 law and his interpolation formula >dilch here fits the 
measiorements very well Indeed. A range of f&lr agreement with the 


V 
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law can lae found and, elmllarly, a range of agreement between the 

correaponding law and the measured .correlation curve (cf. fig. 8 

and reference 22),.. .The range and degree of -agreement is quite similar 
to the one found in the meantime by Townsend (reference 26) at large 
Reynolds numbers. However, in view of the inaccuracy of these measure- 
ments and the arbitrariness in fitting these power laws to measured 
points, agreements such as that presented in figures 7. and 8 should 
be accepted with reserve. Figures 9 to 11 show the change in the cor- 
relation curve with Reynolds number. During the past year the accuracy 
of spectrum measurements has been very much improved and it is hoped 
that future measurements ab large Reynolds numbers can be made with more 
definite results. 


Spectrum Measurements in 20- Inch Tunnel 


Measurements of the spectrum of tiurbv^nce in the 20- inch (correl- 
ation) tiinnel were made with a very much improved technique. Figure 12 
first shows a typical spectrum plotted in the conventional logarithmic 
scale to allow a comparison of the amount of scatter with earlier measure 
ments. Then a series of-spectrums plotted in the reduced form of 

a = F(n)/F(0) against ^ = ^(0)n is shown in figures 13 to l6 for 


various speeds U and grid sizes. All these measurements were made in 

the initial phase of decay. In each figure the spectrum . = (i . if 
is shown for comparison. The meas\ired values lie evidently quite close 
to the simple spectrum for valties of ^ less than about 6. For higher 
values the meas-ured values fall below the simple spectrum, as has been 
discussed in the section "General Considerations." In any case the 
largest part of tte turbulent energy is doubtless contained in a region 


where a=(l+^/^ is a very good approximation. The value of F(0) 
had to be found from the experimental data. This was done by fitting 


a curve of the form — . ■ to the measvirementB of the not-normallzed 

B + - 

spectrum. In this way one finds a value e(0) from which F(o) is 
obtained by normalization, that is, by the division of e(0) by the 




area / e(n) dn. The correctness of e(0) could be checked somewhat by 
comparing the area ^ with the result of a measurement of the 


the mean square of the total output (i,e., ^without frequency analysis) 
which should be equal to the area. This was checked and found to be 
true within l3 percent in all cases. 
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Moments of Spectrum 

The spectrum measurements in the 20 -inch tunnel were accurate enoTigh 
to allow the computation of the second and fourth moments of the spectrum 
density F(n), that is, the quantities 


and 


m n%(n) dn 




Ooo 


nS’(n) dn 


The second moment is used to dete rmine as has been mentioned above; 

the fourth moment furnishes 


I 


G = X^f^"^( 0 ) 


A set of curves showing n^FCn) and n^(n) is given in figures 17 
to 22 . It is seen that the accuracy of the meastirements and the fre- 
quency range is sufficient to evalviate the moments by graphical Integra 
tlon. The resuiLts are tabulated below; 


u 

(cm/sec) 

M 

(cm) 

X 

(cm) 

X2 

(cm?) 

G 

C 

1130 

1.68 

I^ 8.5 

0.05^ 

13.2 

2.6 

1130 

1.68 

107.5 

.123 

13.8 

3.0 

1130 

1.68 

160.0 

.172 

11.9 

2.7 

630 

1.68 

49.5 

.088 

10.7 

3.0 

630 

1.68 

109.5 

.229 

10.2 

3.00 

630 

1.68 

156.5 

.321 

9.39 

3.1 

300 

1.68 

108.5 

.1*63 

8 . 1 ^ 


1500 

1.68 

108.5 

.057 

15.6 


1130 

.818 

108.5 

.150 

10.22 


1130 

.lt -22 

108.5 

.133 

8.20 



During decay G is nearly constant in the initial range and then 
decreases. This agrees with the results of Batchelor and Townsend 
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The measured absolute values for 6 here appear slightly higher. In 
the last column in the table above, the quantity 



is included and C is proportional to LoitziansTcii’s invariant in the 


u2/uo® 


case where the simple spectrum applies. Of course, F(0) and 
vary considerably d\iring decay. As seen from the table, C is nearly 
constant, that is, does not- show a definite trend. Clearly, the accuracy 
in the determination of C is not very good. 


Behavior of Spectrxm at High Frequencies 

It“has been mentioned before that Heisenberg's theory of the spectrum 
of isotropic turbulence leads to an n"' law for high frequencies* It is 
therefore of interest to compute the sixth moments of the spectrum 

n^ = n^(n) dn 

If F(n) ~ n”*^, then the sixth moment should diverge. Plots of n^(n) 
for two speeds etre presented in figures 23 and 24. For the lower speed, 

U = 630 centimeters per second, the area appears to be finite; for the 
higher speed, U - 1130 centimeters per second, the frequency range of the 
measurements is evidently not large enough to decide at all whether or 
not the curves approach zero. 

However, it should be kept in mind that the turbulent energy in the 

frequency range that matters in the n^(n) curve is extremely small. 
Hence, the electrical noise and the pickup in this range are of the order 
from 20 percent to^80 percent of the total output and the measuzed values 
of F(n) are therefore the difference between two fairly close nxjmbers. 

This difference is then multiplied by 10^® to 10^^ in the Important range 
and the resulting curve is therefore not very trustworthy. The only 
definite result which can be asserted so far with confidence is that F(n) 

does not vary more slowly than n'^i^ in the high frequency range. 

Batchelor has recently reported during a visit at GALCIT that Townsend 
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has obtained measurements of 



Since 



Is equivalent to n^, the sixth moment from 

the spectrums. Townsend finds the sixth moment finite. This agrees with 
the present measurements of the spectrum within the limitations discussed 
above. Batchelor Inteirrets this result as confirming the n""^ power 
law plus a rapid cut-off at high frequencies. This conclusion appears 
premat\ire hecaiise of the difficulties stressed above. 


a measurement of 



f Spectrum at Very Large Distances from Qrld 

Very far downstream from the grid. In the so-called final stage of 
decay, the spectrum and the correlation curve should approach a Gaussian 
distribution. Measurements of the time correlation curve were made by 
Batchelor and Townsend (reference 17) and the conclusion was reached 
that the final stage of decay Is rapidly approached for x/M > 300. 

The sx>ectrum of turbulence was therefore measured at x/m > 1000 
behind a fine mesh grid. The measurements are plotted In figure 2^ 
compEored with the Gaussian spectrum compatible with the measured value 
of X; that Is, the spectznm was evaluated as usual and X was deter- 
mined. Then the Gaussian spectrum was plotted which would give the 
same X. 


The measured points at very low frequency are, as mentioned before, 
not too trustworthy. The rest of the cmrve approaches the Gaussian. 
Actually this latter fact becomes more obvious if one tries to plot these 
measurements In the same form as the spectrums at lower values of x/m. 
Figure 26 shows that type of plot and evidently little similarity with 

the spectrum a = (l + ^ can be found here. Similarly, one can 

convince oneself that with Increasing x/M the measured spectrums tend 
toward the Gaiisslan. Still, from the present set of measurements the 
Gaussian form Is not refiched within the accuracy of measurement. Here 
future work Is needed to see whether x/m ■ 1000 Is still too small or 
whether other factors enter. 
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Spectrum of Derivative 

The function F(n) 1b defined as the power spectrum of u(t) . 
Similarly, one can define a power spectrum W(n) of the derivative 
8u/8t. If W(n) is not normalized one has 




W(n) dn 


On the other hand, 


j ’ . g poo 

j * 4it^ n^(n) dn 

u2\8ty Jo 


and, consequently, 

WM - 4jt2n%(n) 


Figure 27 shows measurfflnents of- W(n) plotted In ccanparlson with n%(n). 
The agreement is seen to he excellent. This is a welcome check on the 
consistency of the measurements. 


Zero Counts 

p 

A sample set of 1 values obtained from the counting of zeros 
using an electronic gate clrcvdt is Included in figures 4 and 5- The 

value of was obtained frcm the number of zeros assuming a Ckiusslan 
distribution for both the function u(t) and the derivative 8u/8t. 

Then 1 ? Is given simply by 



Earlier measurements \islng the photocell method gave essentially the 
same result. As computed from the zeros, ^ _ behaves as functions 
of X quite similar to the values obtained from the other methods, 
but the absolute values of X^ as obtained from the counts are consist- 
ently larger, especially at the hl^ speed. How much of-this difference 
in absolute value can be traced to the deviation from the Gaussian 
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distribution and Independent probabilities of u(t) and Su/5t remains 
to be Investigated. That Ng Is nearly Gaussian In | Is shown In 

figure 28. This Is, however, obviously only necessary but not sufficient 
for Independence between u(t) and ^/8t. It was first believed that 
the better agreement for larger values of X reflected the lack of 
resolution of the system. However, the resolution of the new circuit 
should be very much better and the difference between the X measure- 
ments and the zero counts probably has to be traced to other reasons. 


CONCLUDING REMAEKS 


The results of an investigation of the spectrum of isotropic turbu- 
lence led to the following conclusions: 

1. The spectirm and correlation functions In the turbulent flow 
behind a grid at large Reynolds numbers (rN % 10^ based on the grid mesh^ 

show a range of frequencies for which the n"5/3 or r"^^^ laws, respec 
tlvely, apply. Karman has shown that these measurements can be repre- 
sented well by an interpolation formula for the spectrum of the form 


where 


r(o) '' ’ 


5 = |F(0)n 


2. The measurements at smaller Reynolds numbers of the order of 10^ 
show that in the Initial stage of decay the bulk of the turbulent energy 
lies In a frequency range In which the spectrum Is closely approximated 
by the simple relation 


lisl . 

r(o) 


(1 - 


corresponding to an exponential correlation function. 

3. The mean square of the frequency ^ or the corresponding micro- 
scale X was measured by differentiation and from the spectrum. The 
measurements agree closely and the slope of X2(x) Is nearly linear In 
the Initial stage of decay. The value of the slope was found to be a 
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little less than lOv/U. The turhulent energy was found to vary down- 
stream according to the law (u’)^ ** l/t for distances up to about 

100 mesh. The quantity G ■ was obtained from the spectrums 

and G was found nearly constant In the Initial stage of decay. It 
decreases in the later stage} G Increases with speed. These results 
are In general agreement~wlth Batchelor and Townsend's measurements. 

h. The sixth moments of the spectrums were con^nited. Here the limit 
of accuracy of the measurements is reached and conclusions have to be 
tentative. The sixth nonents appear to be finite and hence F(n) at 
very high frequencies should vary ftister than n"T. Certainly F(n) 
does not vary slower than n"*^. 

5. The spectrum at x/M => 1000 was found to approach a Gaussian 
curve. 


6. Using some simplifying asstnnptlonsj X was also obtained from a 
count of the average number of zeros of the fluctuating velocity compo- 
nent. The general trend of X^(x) agrees with the measurements by other 
methods} the absolute values of X^ as obtained from the zero counts are t 
somewhat larger. It was not possible yet to decide whether this differ- 
ence is significant. 


California Institute of Technology 

Pasadena ; Calif., August l6, 19^9 





NACA TN 2l^73 


31 


heperences 


1. Taylor, G. I.: Statistical Theory of Turbulence. Parts I-IV. Proc. 

Roy. Soc. (London), ser. A, vol. 151, no. 873, Sept. 2, 1935, 
pp. 421-478. 

2. Von E^rmlln, Th. : The Fundamentals of the Statistical Theory of 

Turbulence. Jour.' Aero. Scl., vol. 4, no. 4, Feb. 1937, PP* 131-138. 

3. De JOirmin, Theodore, and Hovarth, Leslie: On the Statistical Theory 

of Isotropic Turbulence. Proc. Roy. Soc. (London), ser. A., vol. l64, 
no. 917, Jan. 1938, pp. 192-215. 

4. Kolmogoroff, A. N. : Dissipation of Energy In the Locally Isotropic 

Turbulence. Con^). rend., acad. scl. DESS, vol. 30, no. 4, Feb. 10, 
1941, pp. 301-305j vol. 31, no. 6 , 1^41, pp. 538-540; vol. 32, 
no. 1, July 10, 194 i, pp. l 6 -l 8 . 

5. Batchelor, G. K. : Kolmogoroff * s Theory of Locally Isotropic Turbu- 

lence. Proc. Cambridge Phil. Soc., vol. 43, pt. 4, Oct. 1947, 
pp. 533-559. 

6 . Onsager, Lars: The Distribution of Energy In Turbulence. Abstract, 

Phys. Rev., vol. 68 , nos. 11 and 12, second ser., Dec. 1 and 15, 

1945 , p. 286 . 

7 . Von Welzsdcker, C. F.: Das Spelrtrum der Turbulenz bel grossen 

Reynoldsschen Zahlen. Zeltschr. Phys., Bd. 124, Eieft 12, Sept. 6, 

1948, pp. 614 - 627 . 

8. Heisenberg, W.: Zur statlstlchen Theorle der Turbulenz. Zeltschr. 

Phys., Bd. 124, Heft 7/l2, 1948, pp. 628 - 657 . 

9 . Taylor, Goeffrey Ingram: The Spectrum of Turbulence. Proc. Roy. Soc. 

(London) ser. A, vol. l64, Feb. 18 , 1938, pp. 476-490. 

10. Batchelor, G. K. : Recent Developments In Turbulence Research. 

Introduction. Proc. Seventh Int. Cong. Appl. Mech. (Sept. 1948, 
London), 1948, pp. 27-56. 

11. Townsend, A. A.: Measiirement of Double and Triple Correlation 

Derivatives In Isotropic Turbulence, Proc. CamObrldge Phil. Soc., 
vol. 43 , pt. 4, Oct. 1947 , pp. 560 - 570 . 

12. Batchelor, G. K. : The Role of Big Eddies In Homogeneous Turbulence. 

Proc. Roy. Soc. (London), ser. A, vol. 195, no. 1043, Feb. 3, 1949, 

pp. 513 - 532 . 



32 


NACA TN 2473 


13. Kamp^. de Ferlet, J.: Le tenaeur spectral de la turbulence bomo^ene 

non Isotrope dans un flulde Incoanpresslble. Coiop. rend., acad. 
scl. (Paris), t, 227, Oct. I8, 1948, pp. 78O-761. 

14. Lin, C . C . : On the Law of Decay and the Spectrum of Isotropic 

Turbulence. Vol. 2. pt. I. Proc. Seventh Int. Cong. Appl. Mech. 
(Sept. 1948, London), 1948, pp. 127-139* 

15* Von KR- rman , Th. , and Lin, C. C.: On the Statistical Theory of 

Isotropic Turbulence. Vol. II of Advances In Applied Mechanics, 

R. von Mlses Th. von Karm^, eds.. Academic Press, Inc. 

(New York), 1951> EP*.1-19.* 

16. Loitslanskll, L. G.: Some Basic Laws of Isotropic Turbulent Flow. 

CAHI Rep. No. 440, 1939. (Available in English translation as 
NACA TM 1079, 1945.) 

17. Batchelor, G. K., and Townaend, A. A.: Decay of Turbulence in the 

Final Period. Proc. Roy. Soc. (London), ser. A, vol. 194, 
no. 1039, Nov. 9 , 1948, pp. 527-543* _ 

18. Von Karman, Th.: Progress in the Statistical Theory of Turbulence. 

Proc. Nat. Acad. Scl., vol. 34, no. 11, Nov. 1948, pp. 530-539* 

19. Simmons, L. F. G., and Salter, C.: Eaiperimental Investigation and 

Analysis of the Velocity Variations in Turbulent Flow. Proc. Roy. 
Soc. (London), ser. A, vol. l45, no. 854, June 2, 1934, 
pp. 212-234. 

20. Wang, M. C., and Uhlenbeck, G. E.: On the Theory of the Brownian 

Motion II. Rev. Modem Phys., vol. 17, nos. 2 and 3, April and July 

1945, pp. 323-342. 

21. Dryden, Hugh L.: Turbulence Investigation at the National Bureau of 

Standards. Proc. Fifth Int. Cong. Appl. Mech. (Sept. 1938, 

Cambridge, Mass.), John Wiley & Sons, Inc., 1939, PP* 362-368. 

22. Rice, S. 0.: Mathematical Analysis of Random Noise . The Bell System 

Tech. Jour., vol. 23, no. 3, July 1944, pp. 2^-332} vol. 24, 
no. 1, Jan. 1945, PP* 46-108. 7- 

23. Batchelor, G. K., and Townsend, A. A.: Decay of Isotropic Turbulence 

in the Initial Period. Proc. Roy. Soc. (London), ser. A, vol. 193, 
no. 1035, July 21, 1948, pp. 539-558. 

24. Frenklel, F. N.: On Third-Order Correlation and Vortlcity in Isotropic 

Turbulence. Quart. Appl. Math., vol. VI, no. 1, April 1948, 

pp. 86-90. 



5 


NACA TN 2473 


33 


25. Batclielor, G. K. , and Townsend, A* A.: A Conment on F. N. Frenklel's 

Note "On Third-Order Correlation and Vorticity in Isotropic Turhu- 
lence." Quart. Appl. Math., vol. VII, no. 1 , April 1949 # P* 120 . 

26. Townsend, A. A.: Ejqjerlmental Evidence for the Theory of Local 

Isotropy. Proc. Cambridge Phil. Soc., vol. 44, pt. 4, Oct. 1948, 
pp. 560-565. 


ft 




I 



^ » A * t. y 



KACA TH' 2^73 


































0 


y’ 










✓ 

y 











ff' 










Q 










/ 









1 

1 

m 

a 










✓ 











■r 




10 20 30 40 


50 60 

X 

M 


OJ 

CT\ 


70 80 00 100 110 120 


Figure 3.- Energy decay 'beOilnd grid. M » 1.68 centlmeterB; 
U «= 1130 centimeters per second. 
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Figure 7»- Energy spectrum. RN ='300^000j 4- Inch grid 
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Figure 12.- Energy spectrum of turbulence produced behind grid. 

U = 1130 centimeters per second; x/U - 8 I; M = 1.27 centimeters; 

F(0) = 4.0 X 10“3 second. 
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Figure 13 •- Reduced epectrufli. U = II 30 centimetera per second; 
M = 1.68 centimeters; RN = 12,1|00. 
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Figure l4.- Reduced spectrum. -U = 63 O centimeters per second; 
M = 1.68 centimeters; RN = 69 OO. 
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Figure 15 .- Reduced spectrum. x/M = 65 ; M = 1,68 centimeters. 
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Figure l6.- Reduced spectrum. U = II30 centimeters per second; 

X = 108 centimeters. 
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Figure 18.- Second momentB. U = 630 centlmeterB per second; 

M = 1.68 centimeters. 
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Figure 19<’- Second Tnomen'ta. x/Hg approzliDa'bely 85. 
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Flgui^ 23.- Sixth moineiitB. U =» II30 centimeters per second; 

M = 1.68 centimeters. 
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PlgTore 2k.- Sixth moments. U = 630 centimeters per second; 

M » 1,68 centimeters. 


Figure 25.- Plot of measurementB of spectrum of turbulence at 

x/M =• 1000 behind a fine r»Bh grid and comparison with e"i 
U =-630 centimeters per second; M = 0.l4l- centimeters. 
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Figure 26.- Plot of measurements of spectrum of turljulence at x/M = 1000 
beMnd a fine mesh grid and coD 5 >rarison with — 1 — . u = 630 centi- 


meters per second; M = 1.4l centimeters. 
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Figure 27.- CooQ>arl 8 on of measurements of W(n) and n%(n) . 

U = 1X30 centimeter B per aecom; x/m = 1360* 
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Figure aB.“ Probability distribution of 
U = 1133 centimeters per i 






